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Abstract 

The equation of state of de-confined quark matter within the MIT bag model is 
calculated. This equation of state is used to compute the structure of a neutron star 
with quark core. It is found that the limiting mass of the neutron star is affected 
considerably by this modification of the equation of state. Calculations are carried 
out for different choices of the bag constant. 
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1 Introduction 

Neutron stars (NS) are among the densest of massive objects in the universe. They are 
ideal astrophysical laboratories for testing theories of dense matter physics and provide 
connections among nuclear physics, particle physics and astrophysics. The maximum 
mass of a neutron star is a subject that several theoretical astrophysicists have tried to 
compute it. Below a certain maximum mass, degeneracy pressure prevents an object col- 
lapse into a black hole. To calculate the maximum mass, we require enough information 
about the composition of the star. Different compositions lead to different equations of 
state (EOS). When nuclear matter is compressed to densities so high that the nucleon 
cores substantially overlap, one expects the nucleons to merge and undergo a phase tran- 
sition to de-confined quark matter. Such a system could be realized in two possible ways: 
(a) complete strange quark matter stars (b) neutron stars with a core of quark matter. 
Glendenning [1] has shown that a proper construction of the nucleon-quark phase transi- 
tion inside neutron stars implies the coexistence of nucleon matter and quark matter over 
a finite range of pressure. This has the effect that a core, or a spherical shell, of a mixed 
quark-nucleon phase can exist inside neutron stars. 

In this work, we calculate the structure of the neutron star with a quark core and also 
strange star, and compare our results with our pervious works in which we investigated 
the NS structure without a quark core [2]. 

2 Quark Matter Equation of State 

For the de-confined quark phase, within the MIT bag model [3], the total energy density 
is the sum of a non-perturbative energy shift B (the bag constant) and the kinetic energy 
for non-interacting massive quarks of flavors / with mass m/ and Fermi momentum 
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kp^ = (7r 2 pj)3 where pf is the quarks density of flavor f: 




x 



Vx 2 + l(2x 2 + 1) - sinh -1 x] + -^(n 2 p)l + B, 



(1) 



where x = hkp/mc , p s = pa = p u = p, p is baryon density, and Eq = E/V. We 
assume in this work that u and d quarks are massless and the s quark has a mass equal 
to m — 150MeV. The bag constant B, can be interpreted as the difference between the 
energy densities of the non-interacting quarks and interacting ones, which has a constant 
value such as B = 55 and 90MeV in the initial model of MIT. Inclusion of perturbative 
interaction among quarks introduces additional terms in the thermodynamic potential 
[4] . We try to determine a range of possible values for B by using the experimental data 
obtained at the CERN SPS [5]. 

According to the analysis of those experiments, the quark-hadron transition takes 
place at about seven times normal nuclear matter energy density ( e = 156MeV fm~ 3 ). 
We assume a density dependent B. In the literature there are attempts to understand the 
density dependence of B [6, 7, 8]; however, currently the results are highly model depen- 
dent and no definite picture has come out yet. Therefore, we attempt to provide effective 
parameterizations for this density dependence. Our parameterizations are constructed in 
such a way that at asymptotic densities B has some finite value _Boo : 



The parameter B = B(p = 0) has constant which is assumed to be B = 400 in this work, 
and P is numerical parameter usually equal to p ~ 0.17/m~ 3 , the normal nuclear matter 
density. B^ depends only on the free parameter B . In order to fix B^, we proceed in 
the following way: 

Firstly, we use the equation of state (EOS) of asymmetric hadronic matter character- 
ized by a proton fraction x p = 0.4 and the UVu + TNI potential. By assuming that the 



B(p) = B oo + (B -B oo )ex P -f3{^) 2 



(2) 
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hadron-quark transition takes place at the energy density e = ll00MeVfm~ 3 , we find 
that hadronic matter baryon density is p t = 0.98/m~ 3 (transition density) and at values 
lower than it the quark matter energy density is higher than that of nuclear matter, while 
with increasing baryon density the two energy densities become equal at this density and 
after that the nuclear matter energy density remains always higher. Eq. (1) for quark 
matter with two flavors u and d reduces to: 



3hc 2 



4 4 
Pi +Pd 



+ B, (3) 



where p d = 2p u = 2p. 

Secondly we determine = 8.99 by putting quark energy density and hadronic 
energy density equal to each other (eq = e\ p=Pt ). 

Finally, we calculate the EOS for the three flavors quark matter using 

3 Mixed Phase 

The hadron-quark phase transition takes place within a range of baryon density values. 
In other words, the fraction of space occupied by quark matter smoothly increases from 
zero to unity when eventually the last nucleons dissolve into quarks. In this case, we have 
a mixture of the hadron, quark and electron background in the system. Glendennings 
construction [1] describes a global division of the baryon number between the two phases. 
The equilibrium conditions, in the case where the geometry of droplets is neglected, are 
those for bulk systems. The neutron star matter is assumed to be stable and charge 
neutral. Thermal effects are not expected to play any important role in neutron star cores. 
We neglect such effects and put the temperature T = 0. The equilibrium conditions for the 
quark matter droplet to coexist with the nucleon medium are that pressure and chemical 



potentials in both phases coincide. We choose pressure as an independent variable. The 
coexistence requires that (Gibbs conditions): 

li n N (P)=li%(P), (5) 

and 

p P (p)=p q P (P), (6) 

where p% and p q N are the neutron chemical potential in the nucleon phase (NP) and the 
quark phase (QP), respectively. Similarly, p P and p q P are the proton chemical potential 
in the respective phases. The strange quark and lepton chemical potentials are dictated 
by the conditions of weak equilibrium 

= Vd, (7) 

and 

fid ~ = (8) 

where p\ is lepton chemical potential. 

Using the semi-empirical mass formula, the energy per particle of nuclear matter can 
be expressed as 

E(p,x) = T(p,x) + V (p) + (1 - 2x) 2 V 2 (p), (9) 
where x = p p /p is proton fraction. The kinetic energy contribution is 

T(p, x ) = l^ {^ 2 P) § [(1 - a:)f + ] • (10) 

The functions V and V 2 represent the interaction energy contributions and we can deter- 
mine them from the results of symmetric nuclear matter (x — |) and pure neutron matter 
(x = 0) [9]. Using our results for the nuclear matter with the UVu + TNI potential, we 
have obtained the following fit for Vq and U 2 : 

V (p) = -559.9p 5 + 1695.62312p 4 - 1946.86437p 3 + 1327.04p 2 - 411.57428p - 0.30327,(11) 



V 2 (p) = 191.21p 5 - 626.712p 4 + 776.623p 3 - 473.40909p 2 + 141.8709p + 4.75638- (12) 



From Eqs. (9) and (10), we obtain the chemical potentials of neutrons and protons as: 

h 2 



2m 



(3tt 2 p) 5 [(l-x)l +x(l-x^]+V (p)+ P Vo(p) 



+ (1 - 2x) 2 pV 2 \p) + (1 - Ax 2 )V 2 {p) + mc 2 , 



(13) 



and 



x°* + (1 - x)x* + V (p) + pVq(p) 



+ (1 - 2x) 2 pV 2 \p) + (-3 + 8x - Ax 2 )V 2 (p) + mc 2 



(14) 



The quark chemical potential with favor / is 



H = 



m 2 f c 4 + h 2 c 2 (n 2 p f )i 



(15) 



From Eqs. (1), (4), (7) and (15), we obtain the chemical potential of quark matter: 



4ttW 



„ „ „ f dD as/ 



where D = |^ [aVx 2 + l(2x 2 + 1) - sinh -1 



We also have 



(16) 



P^r — 2p d + 



(17) 



and 



p% = 2p u + p d - 



(18) 



By plotting p P versus p^ for both nucleon and quark phases, we can identify the cross 
point of two curves that satisfy the Gibbs condition. As the chemical potentials determine 
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the charge densities of the two phases, the volume fraction occupied by quark matter, x, 
can be obtained by exploiting the requirement of global charge neutrality: 

XP c q + (l-x)Pp~Pe = 0, (19) 

where p c is the quark charge density. The total energy density and density of mixed phase 
(MP) are given by: 

£mp = X£qp + (1 - x)£np, (20) 

and 

PMP = XPQP + (1 - X)PNP- (21) 

We plot pressure versus baryon density for hadron, mixed and quark phases in Figures 1 
and 2, for bag constants B = 90 and density dependent B, respectively. It is seen that 
there is a mixed phase at a range of densities. A pure quark phase is presented at densities 
beyond this range and a pure hadronic phase is presented at densities below it. 



4 Structure of Neutron Star with a Quark Core 

We calculate the structure of a neutron star for various values of the central mass density, 
e c , by using the equation of state and numerically integrating the general relativistic 
equation of hydrostatic equilibrium, Tolman-Oppenheimer-Volkoff (TOV) equation[10]. 
The derivation of TOV-equation can be found in standard textbooks [11, 12, 13, 14]. 
For a neutron star with a quark core, we use the following equations of state: 

• Below the density 0.05/m -3 , we use the equation of state calculated by Baym et al. 
[15]. 

• From this density up to the beginning point of the mixed phase, we use the equation 
of state which is calculated with the UVu + TNI potential [2]. 
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• For the range of densities in the mixed phase, we use the equation of state which 
was calculated in the previous section. 

• Beyond the end point of the mixed phase, we use the equation of state of pure quark 
matter which was calculated in section 2. 

Calculations are done both for constant B = 90, and density dependent B. Based on 
these EOS's, we calculate the mass and radius of the NS with quark core. The calculations 
are also repeated for the strange star (i.e. pure quark matter). We plot the NS mass 
versus central mass energy density for B=90 and density-dependent B, in Figures 3 and 
4, respectively. The NS mass versus radius for quark core NS and strange star are plotted 
for B=90 and density-dependent B in Figures 5 and 6, respectively. For the sake of 
comparison, we have also plotted our previous results of the neutron star structure without 
quark core, in these figures. It is seen that there is a profound difference between the new 
results for NS with a quark core and those of NS without a quark core. 

The extracted maximum mass of a NS and the corresponding radius and central mass 
density for both cases B=90 and density dependent B are presented in Tables 1 and 
2, respectively. It is seen that the inclusion of the quark core leads to a considerable 
reduction of the maximum mass, while the radius is not affected appreciably. Note that 
the maximum mass for the NS with quark core and B=90 is quite near to the observed 
maximum mass of neutron stars [16]. 

The maximum mass energy density versus the radial coordinate for NS without core, 
NS with a quark core and strange star are plotted in Figures 7 and 8 for B=90 and density- 
dependent B, respectively. It can be seen that a major part of the core is composed of 
pure quark matter (about 8 Km). A layer of mixed phase (thickness about 1.5 Km) exists 
between the core and a thin crust. 
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5 Summary 



As we go from the center toward the surface of a neutron star, the state of baryonic matter 
changes from the de-confined quark-gluon to a mixed state of quark matter and hadronic 
matter, and thin crust of hadronic matter. The transition between these states occurs in 
a smooth way. In order to calculate the structure and the mass limit of neutron stars, it 
is important to have a fairly accurate physical description of these states. 

In this paper, we calculated the equation of state of the de-confined quark phase within 
the MIT bag model. We then calculated the mixed phase of nucleons and quarks. The 
equilibrium volume fractions of nucleon and quark matter were obtained by applying the 
Gibbs condition. Curves were presented which showed the dependence of pressure on the 
baryon density. 

Our results for the equation of state were then used to calculate the structure of a 
neutron star with a quark core. As usual, the Tolman-Oppenheimer-Volkoff equation were 
integrated from the center to the surface of the neutron star where the density drops to 
zero. Calculations were carried out both for B=90 and a density-dependent B. 

The maximum mass, radius, and central mass density of neutron stars with a quark 
core and strange star were calculated and compared with the traditional neutron star. It 
was found that the limiting mass decreases when the quark core is taken into account. 
This brings the maximum mass closer to the observational limits. 
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Table 1: Maximum gravitational mass (M max ), corresponding radius(R) and central mass 
density(e c ) for B = 90. 



star 


M max (M ) 


R(Km) 


e c (10 u gr/cm 3 ) 


NS 


1.98 


9.81 


27.17 


NS+quark core 


1.57 


9.73 


33.27 


stange star 


1.34 


7.77 


34.81 



Table 2: Maximum gravitational mass (M max ), corresponding radius(R)and central mass 
density(e c ) for density dependent B. 



star 


M max (M ) 


R(Km) 


e c (10 u gr/cm 3 ) 


NS 


1.98 


9.81 


27.17 


NS+quark core 


1.75 


9.66 


28.92 


stange star 


1.63 


8.2 


28.92 
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Figure 1: The pressure versus baryon density for Hadron Phase (solid line), Mixed Phase 
(dashed line) and Quark Phase (dotted line) for B=90. 
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Figure 2: The pressure versus baryon density for Hadron Phase (solid line), Mixed Phase 
(dashed line) and Quark Phase (dotted line) for density dependent B. 
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Figure 3: The gravitational mass versus central mass density for different cases with bag 
constant B=90. 
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Figure 4: The gravitational mass versus central mass density for different cases with 
density dependent B. 
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Figure 5: The mass-radius relation for different cases with bag constant B=90. 




Figure 6: The mass-radius relation for different cases with density dependent B. 
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Figure 7: Mass density as a function of radial coordinate for neutron star (dotted line) 
neutron star with quark core (solid curve) and strange star (dashed curve) with B=90. 
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Figure 8: Mass density as a function of radial coordinate for neutron star (dotted line) 
neutron star with quark core (solid curve) and strange star (dashed curve)with density 
dependent B. 
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